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Abstract 

In this paper Fokker-Planck-Kolmogorov type equations associated with stochastic 
O ,1 differential equations driven by a time-changed fractional Brownian motion are derived. 

Two equivalent forms are suggested. The time-change process considered is either the 
first hitting time process for a stable subordinator or a mixture of stable subordinators. 
A family of operators arising in the representation of the Fokker-Plank-Kolmogorov 
equations is shown to have the semigroup property. 



^ 1 Introduction 

On 

In this paper we establish Fokker-Planck-Kolmogorov type equations associated with stochas- 
tic differential equations driven by a time-changed fractional Brownian motion. A (one- 
dimensional) fractional Brownian motion Bf 1 is a zero-mean Gaussian process with contin- 
uous paths and correlation coefficient 

- C B {s, t) = E(B^Bf) = \{s™ + e H -\s- tn, (1) 

h : 

where the Hurst parameter H takes values in (0,1). If H = 1/2, then the correlation 
disappears, and Bf becomes a standard Brownian motion. Stochastic processes driven by 
a fBM are of increasing interest for both theorists and applied researchers due to their wide 
application in fields such as mathematical finance, astrophysics, turbulence, etc. 

Fractional Brownian motion Bf, like standard Brownian motion, has nowhere differen- 
tiable sample-paths. The covariance between increments over non-overlapping intervals is 
positive, if | < H < 1, and negative, if < H < |. Increments of Bf 1 exhibit long range 
dependence if | < H < 1. The Hurst parameter H can be extended to H = 1 as well, the 
corresponding fBM having the form B\ = tN, where N is the standard normal random vari- 
able. If H e (0, 1/2) U (1/2, 1), then B? has the representation pj[25] B? = jj K H (t, s)dB s , 
where Kn{t,s) is expressed through Gauss' hypergeometric function. FBM is not a semi- 
martingale unless H = 1/2, so the usual Ito's stochastic calculus is not valid. Nevertheless, 
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there are several approaches [31 El El 125] to a stochastic calculus in order to interpret in a 
meaningful way a SDE driven by an m- dimensional fBM Bjr of the form 

X t = X + [ b(X s )ds+ [ a{X s )dB*, (2) 
Jo Jo 

where mappings b : MJ 1 — > IR n and o : K n — > MJ lxm are Lipschitz continuous and bounded; X 
is a random variable independent of Bf . We do not discuss in this paper these approaches 
referring the interested reader to El 125]- Instead, we focus our attention on the FPK 
equation associated with SDE ((2|) driven by fBM whose form is [2j [9] 

^M^, H ^ + ^E, W ^ (3) 

3=1 3 j,k=l 3 

with the right hand side dependent on the time variable t, which, in fact, reflects the presence 
of correlation. Functions a,jk(x),j, k = 1, n are entries of the matrix A(x) — a(x) x <r T (x), 
where a T (x) is the transpose of matrix a{x). By definition A(x) is positive definite: for 
any and (6l" one has YTj k=i a jk{x)^k > C|^| 2 , where C is a positive constant. 

Additionally, u(t, x) in equation (E]) satisfies the initial condition 

u(0,x) = ip(x), x e W 1 , (4) 

where <p(x) belongs to some function space, or is a generalized function. In the particular 
case of FPK equation associated with SDE (J2J), <f(x) = fx (x), the density function of X . 
If X = xq G W 1 , then f(x) = 5 X (x), Dirac's delta with mass on xq. In this case the solution 
to the FPK equation is understood in the weak sense. 

In the sequel we use the following conventional notation: 

- d - d 2 

B{x,D x ) = ^b J {x) — , A{x,D x )=^a jk (x)-Q-Q-, (5) 

j=l 3 j,k=l 3 

and 

L 7 {t, x, D x ) = B(x, D x ) + 2±l t J A (x, D x ), (6) 

where 7 = 2H — 1. Due to the condition on the coefficients a^ix) stated above, A(x, D x ) is 
an elliptic operator. If 7 = 0, or equivalently H = 1/2, then the operator L (t,x, D x ) takes 
the form (coefficients not depending on t) 

L {t, x, D x ) = L{x, D x ) = B(x, D x ) + ^A{x, D x ), (7) 

and equation (J3j) coincides with the FPK equation associated with the SDE driven by 
Brownian motion B t (see, e.g. [29J) 

^^- = L(x,D x )u(t,x). (8) 
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The fractional FPK equation is obtained from equation (jSj) upon replacing the first order 
derivative on its left hand side by the time-fractional derivative D% in the sense of Caputo- 
Djerbashian [10]. By definition, the Caputo-Djerbashian derivative of order /3 is given by 

Df ^ = f(i^yjf^' ^ <1 ' (9) 

where T(-) stands for Euler's gamma function. Introducing the fractional integration operator 

J a f(t) = -i- fit - r) a -V(r)rfr, a > 0, 

r(«) J 

one can represent D@ in the form D% = J 1- ^^- We also write D^ t emphasizing that the 
fractional derivative acts with respect to the variable t. An equivalent but slightly different 
representation of the fractional FPK equation is possible through the Riemann-Liouville 
derivative also, see e.g. [28]. The obtained Cauchy problem for the time- fractional FPK 
equation 

D%v(t, x) = L(x, D x )v(t, x),t>0, x e R n , (10) 
v{0,x) = <p(x), x G M n , (11) 

describes the dynamics of a stochastic process driven by a time-changed Brownian motion 
(see Section [2]). Such equations appear in many fields, including statistical physics [2H 153] , 
finance [13], hydrology [I], cell biology [25], etc. Existence and uniqueness theorems related to 
the Cauchy problem for fractional differential equations, as well as more general distributed 
order equations, can be found in [HI [161 ED] . Instead, we focus on how fractional order 
FPK equations are obtained from non-fractional FPK equations. 

By definition, a time-change process is a stochastic process with continuous nondecreasing 
sample paths starting at 0. Let B be a standard Brownian motion and E be the time-change 
process given by the first hitting time process for an independent stable subordinator with 
index /3. If one replaces the driving process B by a composition B o E, then the left hand 
side of equation (jSj) becomes the fractional derivative of order /3, and the right hand side 
remains unchanged. For details we refer the reader to [HI |TH]. As we will see, this is not 
the case for fractional FPK equations associated with SDEs driven by time-changed fBM 
(Section [3]). 

Fractional FPK equations associated with SDEs driven by a time-changed fBM (see 
equation (fT2"j) below) have not yet been determined. Meerschaert et. al. [22] studied the 
continuous time random walk (CTRW) limits for certain correlated random variables, which 
include linear fractional Levy stable motions, and in particular, fractional Brownian motion. 
For the latter, the scaling limits represent time-changed fBM, where the time-change process 
is the inverse to a stable subordinator. Authors of that paper write, 'An interesting open 
question is to establish the governing equation for the CTRW scaling limit." A particular 
case of our Theorem 3.1 answers that question. 

There are several approaches for deriving equation (fTUj) . including via semigroup theory 
[Tj, [H], master equations [201 EZ], and continuous time random walks [TTj [T21 [231 E2]- In 
this paper we use a different technique, which can be extended for equations with a time 
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dependent right hand side as well, including equations of the form flS]). This technique is 
close to the method used in [T7] . 

SDEs driven by fBM are studied by several authors using different approaches; for refer- 
ences we refer the reader to [5]. SDEs driven by time-changed Brownian motion are discussed 
in [T5]. The associated fractional FPK equations driven by time-changed Levy processes 
when the time-change process is the inverse to an arbitrary mixture of stable subordinators 
are studied by Hahn, et. al. in [TJ]. Note that any time-changed semimartingale is again 
a semimartingale. However, since fractional BM is not a semimartingale if H ^ 1/2, the 
methods used in [TJ] and [15] are not applicable in this case. We plan to discuss a possible 
interpretation of SDEs driven by a time-changed fractional Brownian and linear fractional 
stable motion in a separate paper. Thus, in the present paper we derive FPK type equations 
associated with the SDE 



where E t is the inverse to an arbitrary mixture of stable subordinators with indices in (0, 1). 
Throughout the paper we assume that E t is independent of the driving process Bf. An 
important particular case is when E t is the inverse to a single stable subordinator. The main 
ideas used in this paper will be illustrated in this simpler case. The associated FPK equation 
can be represented as a time-fractional order differential equation, but the right hand side 
does not coincide with the right hand side of equation (|2D, unless 7 = (or equivalently, 
H = 1/2). However, in the case of zero drift (i.e. b(x) = 0), the FPK equation can be 
obtained with the same operator as on the right hand side of ([3]), but in this case the left 
hand side is not a time-fractional differential operator. This difference of FPK equations is 
an essential consequence of the correlation of the increments of the fBM that is the driving 
process of the corresponding SDEs. 

The paper is organized as follows. Section 2 illustrates the method of this paper when 
the driving process is a time-changed Brownian motion. The results obtained in this section 
further clarify properties of density functions of processes which are inverses of arbitrary 
mixtures of stable subordinators. In Section 3, two equivalent FPK equations associated 
with SDEs driven by time-changed fBM are obtained extending the technique used in Section 
2. Furthermore, the family of operators appearing in the FPK equations is shown to have 
the semigroup property. 

2 FPK equations associated with SDEs driven by a 
time-changed Brownian motion 

Consider a SDE driven by a time-changed Brownian motion 



where b(x) and a(x) are Lipschitz continuous mappings and E t is the first hitting time process 
for a stable subordinator W t with stability index (3 G (0, 1). The process E t is also called an 




(12) 



X t = x + [ b(X s )ds+ [ a(X s )dB Ea , t > 0, 



(13) 



Jo Jo 
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inverse to W t . The relation between E t and W t can be expressed as E t = min{r : W T > t}. 

i 

The process W t , t > 0, is a self-similar Levy process with = 0, that is W ct = cPW t 
as processes in the sense of finite dimensional distributions, and its Laplace transform is 
K(e~ sWt ) = e~ tsP . The density fw 1 ( T ) of W\ is infinitely different iable on (0,oo), with the 
following asymptotics at zero and infinity 150] : 

o 2-/3 

2(1-/3) _ ^ 

fa ^ ~ r(i-V 1+ ^ r °°' (15) 

Since Wt is strictly increasing, its inverse process E t is continuous and nondecreasing, but 
not a Levy process. Likewise the time-changed process B Et is also not a Levy process (see 
details in 1141). The associated FPK equation in this case has the form 

D?v{t,x) = L{x,D x )v{t,x), (16) 

with the initial condition v(0, x) = S x (x), where L(x,D x ) is defined in (J7|), and is the 
fractional derivative in the sense of Caputo-Djerbashian. 

Notice that solutions to equations (1161) and (jBD are connected by a certain relationship. 
Namely, a solution v(t,x) to equation (ITOj) satisfying the initial condition PJ can be repre- 
sented through the solution u(t,x) to equation (jSJ), satisfying the same initial condition (J3J), 
by the formula 

/■oo 

v(t,x)= f t (r)u(r,x)dr, (17) 



Jo 

where /t(r) is the density function of E t for each fixed t > 0. If (i) is the density function 
of Wi, then 

/.w = 4*. <^> = -I i ^ ^. = (i8) 

Since fw (u) G C°°(0, oo), it follows from representation ( fl8|) that / t (r) G C°°(R+), where 
= (0, oo) x (0, oo). Further properties of /t( r ) are represented in the following lemma. 

Lemma 2.1 Let /t(r) 6e the function given in ( fl8l) . T/ien 

(oj lim^ + o/t( r ) — ^o( r ) i/ie sense of the topology of the space of tempered distributions 
V'(R); 

(b) lim T ^ +0 /t(r) = ff^y, t > 0; 

fcj lim^oo f t (r) = 0, t > 0; 

(d) C t ^ s [f t (r)](s) = s^e"™', s > 0, r > 0, 

where Ct->s denotes the Laplace transform with respect to the variable t. 
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Proof, (a) Let ip{r) be an infinitely differential) le function rapidly decreasing at infinity 
We have to show that lim^ +0 < ft^ip >= ^(0). Here < f t ,i/j > denotes the value of 
ft E V(R) on ip. We have 

POO 

lim < / t (r), V(r) > = lim / / t (r)^(r)dr 



lim 



f Wi {u)^({-f)du 

1 It 



= V(0) / /w»du = V(o)- 

Parts (6) and (c) follow from asymptotic relations (fT5|) and (|T4|) . respectively. Part (d) is 
straightforward, just compute the Laplace transform of ftij). ■ 

Due to part (b) of Lemma \2.1\ f t G C°°(0, oo) for each fixed r > 0. Hence, the fractional 
derivative D^ t f t (r) in the variable t is meaningful, and is a generalized function of variable 

T. 



Lemma 2.2 Function /t(r) defined in ({TBI) /or eac/i t > satisfies the following equation 

Mr), (19) 

Ul L \L — p) 

in the sense of tempered distributions. 



d t~ 
D *Jt(r) = -^ft(r)-^ r _ f3y 



Proof. The Laplace transform (in variable t) of the left hand side, using the definition 
( ITS]) of /t(r), equals 



C t ^M,Jt(r)](s) = s' J C t ^[f t (r)}(s) - Hm f t (r) = s^e - - , 



On the other hand, the Laplace transform of the right hand side, 



t— >s 



d t-P 



6 n (r) 



d 2 ,1 



Or 2 s 



-e 



W-i*-"" - sP-Xir), s>0, 



completing the proof. 



Derivation of fractional FPK equation. Now it is easy to show the derivation of the 
fractional order FPK equation (fT6|) . a solution of which is given by v(t, x) in (ITT]) . We have 



/"OO 

D? it v(t,x)= / D% t U 
Jo 



t(r)w(r, x)dr 



u(t, x)dr 



+ 



- lim [ft{r)u(T, x)) + lim[/ t (r)u(r, x)] 

T— »oo t— >0 
oo £-/3 
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Due to Lemma 12.14 P ar t (c) implies the first term vanishes since u(t,x) is bounded, while 
part (b) implies the second and last terms cancel. Taking into account OH]), 



D p t v(t,x)= / f t (r)L(x,D x )u(T,x)dr = L(x,D x )v(t,x). (20) 
J o 

Moreover, by property (a) of Lemma 12.11 

lim v(t,x) —< 8o(t), u(t, x) >=u(0,x) = S x (x). m 

This technique extends to the more general case when the time-change process is the 
first hitting time for an arbitrary mixture of independent stable subordinators. Let p(s) = 

s /3 <i/i(/3), where pi is a finite measure with supppi C (0,1]. Let Wf 1 be a nonnegative 
stochastic process satisfying E(e~ sW/(M ) = e~ tp ^ s \ and E£ = min{r : W£ > t}. The process 
Wf represents a mixture of independent stable subordinators with a mixing measure \x (see 



Theorem 2.3 Letu(t,x) be a solution of the Cauchy problem 

du ^ x ) = L (x, D x )u(t, x), t > 0, iGR", (21) 
u(0,x) = (p(x), x e W 1 . (22) 

Then the function v(t,x) = J °° f£ (t)w(t, x), where ftij) is the density function of Ef, 
satisfies the Cauchy problem 

D^v(t,x) = [ D^ t v(t,x)d/i{^) = L(x,D x )v(t,x), t > 0, x G MP, (23) 
Jo 

v(0,x) = f(x), x e W 1 . (24) 

The proof of this theorem requires two lemmas which generalize Lemmas 12.11 and 12.21 
Define the function ^ 

Mt) = J frrT) d ^ )l t > °' (25) 

Lemma 2.4 Let ft(r) be the function defined in Theorem \2.3[ Then 

(a) lim^+o ft(r) = 5 (r), r > 0; 

(b) lim T ^ + o// 1 (r) = $ M (t), *>0; 

(c) lim T _ >oo // t (r) = 0, t>0; 

(d) C t ^[ft(T)](s) = fe-*), s > 0, r > 0. 
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Proof. First notice that /f(r) = /^(r) — — ~§pJfw?{t)-> where J is the usual integration 
operator. The proofs of parts (a) — (c) are similar to the proofs of parts (a) — (c) of Lemma 
12.11 Further, using the definition of Wf 1 , 

C t ^[ft(r))(s) = -~U+.\Jw?im*) = ^e**>\ s > 0, 
which completes the proof. ■ 

Lemma 2.5 The function //*(t) defined in Theorem \2.'J\ satisfies for each t > the follow- 
ing equation 

D,,tft(r) = -§pf?(?) - S (r)%(t) } (26) 
in the sense of tempered distributions. 

Proof. Integrating both sides of the equation £t_>. s t /f(r)] = s /3 ' Ct^ fS [fti T )}i. s ) ~ 
s^ _1 5 (r), and taking into account part (d) of Lemma EH yields 



C 

It is easy to verify that the latter coincides with the Laplace transform of the right hand 
side of ([26]). ■ 

Proof of Theorem \2.3l Using Lemma 12.51 we have 

POO 

D^vit, x) = / D^ft(r)u(r, x)dr = - lim [ft(r)u(r, x)} + lim[/f (t)u(t, x)\ 

f°° d f°° d 

+ J o fnT)-^u(T,x)dr-%(t)u(0,x) = J ft(r)—u(T,x)dr, 

since all the limit expressions vanish due to parts (b) and (c) of Lemma I2T41 Now taking into 
account equation f[2"Tj) . 

POO 

DiAt,x)= / ft t (T)L(x,D x )u(T,x)dr = L(x,D x )v(t,x). (27) 
Jo 

The initial condition (122]) is also verified by using property (a) of Lemma 12.41 
lim v(t, x) =< So(t), u(t, x) >= u(0, x) = (f(x). m 

3 FPK equations associated with SDEs driven by time- 
changed fBM 

Now let us focus on the FPK equation associated with the SDE (1121) driven by a time-changed 
fBM Bfi t . Recall that the FPK equation associated with a SDE driven by a fBM (without 
time-change) has the form [2j [9] 

^^l = Ly (t,x,D x )u(t,x), (28) 
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where L^(t,x, D x ) is defined in fl^]) and the Hurst parameter H is connected with 7 via 
2H — 1 = 7. Again for simplicity, we consider a time-change process E t inverse to a single 
stable subordinator W t , though mixtures of stable subordinator can be treated similarly. 
Hence, the density function ftij) of E t possesses all the properties mentioned in Lemmas 
O and El 

Theorem 3.1 Let u(t,x) be a solution to the Cauchy problem 

= B {x, D x )u{t, x) + ^±±t^A(x, D x )u{t, x), t > 0, x E W 1 , (29) 
u(0,x) = (p(x), x e R n . (30) 

Let /t(r) be the density function of the process inverse to a stable subordinator of index 
f3. Then v(t,x) = J °° ft(T)u(r,x)dr satisfies the following Cauchy problem for a fractional 
order differential equation 

D?v(t, x) = B(x, D x )v(t, x) + ^-^G 7;t A(x, D x )v{t, x), t > 0, xeW 1 , (31) 
v(0,x) = (p(x), x e R n , (32) 
where the operator G lyt o-cts on the variable t and is defined by 

1 rC+ioo ~f \ 

G^ t v(t,x) = {Jr(7 + l)j}- f, C:\ t [— J c ^ {sP _J )l+l dz]{t), (33) 

with < C < s, and z 13 = ePLn(z) ^ j\, n ^ z ^ being the principal value of the complex \n(z). 
Proof. Using the properties of /t(r) we obtain 

DlMt, x) = B{x, D x )v(t, x) + ^— — A(x, D x ) / / t (r)r^(r, x)dr 



7 + 1 

= B(x,D x )v(t,x) H — A(x,D x )G lit v(t,x), 

where 

poo 

O r>t v{t,x)= / f t {r)^u{r,x)dr. (34) 
J 

It follows from the definition ( TlTl) of v (t, x) that if 7 = 0, then G ,t is the identity operator. 
To show representation fl33|) in the case 7 7^ we find the Laplace transform of G J}t v(t,x). 
In accordance with the property (d) of Lemma ( 12. ip we have 

POO 

C[G^ t v(t, x)](s) = s^ 1 / e~ TS Vu(T, x)dr = s^C^u^, x)](s p ). 

Jo 

Obviously, if 7 = 0, then £[G jV (t, x)](s) = s /3_1 £t(s /3 , x), which implies v(s, x) = s /3_1 £t(s /3 , x). 
If 7 ^ 0, then 

C[tMt,x)}(s) = C[r](s)*u(s,x) = ^- [ C f il+ ^ 1 u(z,x)dz } (35) 

2m Jc-ioc I s - Z ) l+ 
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where * stands for the convolution of Laplace images of two functions, and < c < s. Now 
using the substitution z = with Ln(£) the principal part of the complex function 

ln(C), the right hand side of (I35j) reduces to 

£[r«(«,,)] W = ± f™ J ^±±e-^, x)iC 

' ^ r <^ +1 » (36) 



The last equality uses the relation v((,x) = C /3_1 ^(C^) x )- Further, replacing s by s 13 and 
taking the inverse Laplace transform in (136]) we obtain the desired representation (133]) for 
the operator Cr 7jf . In accordance with part (a) of Lemma [2TT1 we have v(0, x) = u(0, x), which 
completes the proof. ■ 

In the more general case when the time-change process E t is the inverse to Wf 1 , the 
mixture of stable subordinators with the mixing measure p, a representation for the FPK 
equation is given in the following theorem. 

Theorem 3.2 Let u(t,x) be a solution to the Cauchy problem ( 129]) -(130]) . Let //*(r) be the 
density function of the process inverse to Wf 1 . Then v(t,x) = f °° f^{r)u{r 1 x)dT satisfies the 
following Cauchy problem for a fractional order differential equation 

D^v(t, x) = B(x, D x )v(t, x) + 1^G$ tt A{x, D x )v(t, x), t > 0, x e M n , (37) 
v(0,x) = <p(x), x e R n . (38) 
The operator G*i t acts on the variable t and is defined by 

G >( t, X) = (,) . [£<I±i) [2 ( ^ Z l%;lM M, (39) 

where * denotes the usual convolution of two functions, < C < s, p(z) = e l3 ^' n ^dp(f3), 
m^z) = & ^ , and $ M (t) is defined in (EHJ) . 

Proof. The proof is similar to the proof of Theorem 13.11 We only sketch how to ob- 
tain representation (139]) for the operator G^ t v(t, x) = J °° fl L {r)T 1 u{r ) x)dT. The Laplace 
transform of G^ t v(t,x) due to part (d) of Lemma [2.4] is 

= ^-C[t^u(t,x)}(p(s)),s > 0. 

Since £[$ M ](s) = s > 0, we have 

Further, replacing s by p(s) in (135]) . and using the substitution z = p(() = / V in ^d//(/3), 
in the integral on the right side of ( 135]) yields the form (139]) . ■ 

The following theorem represents the general case when the time-change process E t is 
not necessarily the first hitting time process for a stable subordinator or their mixtures. 
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Theorem 3.3 Let 7 G (—1, 1). Let E t be a time-change process and assume that K(t,r) = 
/b (t) satisfies the condition: lim T ^ +0 [(^K^t, r)l < 00 /or a// 1 > 0. Lei ii^ fre an operator 
acting in the variable t such that 



H t K^T) = -^-[K(t,T)( t -y 



6 (r) lim 

r->+0 



f '~VK(t,T) 
T 



(40) 



TTien i/ie function v(t,x) = J °° K(t,r)u(T,x)dr, where u(t,x) is a solution to the Cauchy 
problem (|29|) -( 130|) . satisfies the equation 

H t v{t, x) = f ( G^ t B{x, D x )v(t, x) + ^J-l^A{x, D x )v{t, x), t > 0, r > 0, (41) 

and t/ie initial condition v(0,x) = u(0,x). Here G_ 7jt f(t, x) = / °° K{t,r)r^' y u(T,x)dT. 

Remark 3.4 Obviously, if 7 7^ 0, then H t can not be a fractional derivative in the sense of 
Caputo (or Riemann-Liouville). A representation of H t in cases when E t is the inverse to a 
stable subordinator, is given in Corollary 13. 51 

Proof. We have 



H t v(t,x) 



H t K(t, t)u(t, x)dr 



d_ 

^9t 



I. 



K(t,r)(-y] +5 (t) lim U-)^K(t,r) \u(r,x)dr 

T J t->+0 L T 



-t 7 lim [K(t, t)t~ j u{t, x)) + t 7 lim [K(t, t)t- j u{t, x)) 



(42) 



+ 



T 



) 7 X(t,r) u(0,x) 



Obviously, the first term on the right of (1421) is equal to zero, since for each fixed t > 
function K(t,r) is bounded when r — > 00 and u(r,x) decays at infinity. The sum of the 
second and last terms, which exist by the hypothesis of the theorem, also equals zero. Now 
taking equation fl29|) into account, 



H t v(t, x) = FB(x, D x ) / K(t, r)r" 7 n(r, x)dr + 



7 + 1 



t J A(x,D x )v(t,x). 



Further, since Eq = it follows that 



\imv(t, x) 



5o(r)u(T,x)dT = u(0,x). 



Let n 7 denote the operator of multiplication by t 7 , i.e. Il 7 /i(t) = Ph(t), h G C(0,oo). 
Applying Theorem 13.31 to the case K(t,r) = /t(r) in conjunction with Theorem 13.11 we 
obtain the following corollary. 



Corollary 3.5 Let 7 < and K(t,r) = ft(r), where ftij) is defined in (fl8l) . Then 
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(it) H t = n 7 c_ 7jtJ Df . 

This Corollary yields an equivalent form for FPK equation (T3TT) in the case when E t is 
the inverse to the stable subordinator with index /3 and 7 < : 

H t v{t, x) = t"G^ tt B(x, D x )v{t, x) + 1±lt^A(x, D x )v(t, x), (43) 

with H t as in Corollary 13.51 

Suppose the operator in the drift term B(x, D x ) = 0. Then equation (143]) takes the form 

H t v{t, x) = l±±VA{x, D x )v{t, x). (44) 

Consequently, given a FPK equation associated to an SDE driven by a non-time-changed 
fBM, the FPK equation for the analogous SDE driven by the time-changed fBM cannot be of 
the form: retain the right hand side and change the left hand side to a fractional derivative. 
Moreover, if a fractional derivative is desired on the left hand side in the time-changed case, 
then the right hand side must be a different operator from that in the non-time-changed 
case. 

Notice that FPK equation (143]) is valid for 7 G (0, 1) as well. Indeed, part (ii) of Corollary 
13.51 can be rewritten in the form G ljt — G- 7 t f° r 7 > 0. For 7 < part (ii) of Corollary 13.51 
also implies (G~j) _1 = Gzh t = G 7 j. Now applying operators G_ 7i j and n 7 consecutively to 
both sides of fl3T} we obtain for all 7 G (-1, 1). 

Analogously, the FPK equation obtained in Theorem 13.21 with the mixing measure [i can 
be represented in its equivalent form as 

H?v(t, x) = t^G^ t B(x, D x )v(t, x) + ^^t 7 A(x, D x )v{t, x), t > 0, r > 0, (45) 

where Hf* = U^G^ t D^. We leave verification of the details to the reader. 

The equivalence of equations f l3T|) and f )43|) and the equivalence of equations fl37|) and 
( 145]) are obtained by means of Theorem 13.31 This fact can also be established with the help 
of the semigroup property of the family of operators {G 7 , — 1 < 7 < 1} : 

POO 

G 7 g(t)= / f t {r)^h{r)dT = 7 1 h{t) 1 (46) 
Jo 

where h G C°°(0,oo) is a non-negative function rapidly decreasing at infinity. Denote the 
class of such functions by U. Functions g and h in (j46p are connected through the relation 
g(t) = J °° f t {r)h{r)dT = J-h(t). It follows from the behaviour of /t(r) as a function of t, that 
g G C°°(0, 00), but not necessarily integrable. On the other hand, obviously, operator J 7 is 
bounded, WJ-'hW < \\h\\ in the sup- norm, and one-to-one due to positivity of /t(r). Therefore, 
the inverse J 7-1 : FU — > U exists. Let a tempered distribution H(t, r) with suppH C be 
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such that T l g{t) = J °° H(t,r)g(r)dr. Since /t(r) G TU as a function of t for each r > 0, 
for an arbitrary h £ U one has 

poo poo 

h(t)=f- l Fh(t)= I H{t,s)( f s (r)h(r)dr)ds 
Jo Jo 

poo 

h(r)( / H(t,s)f s (r)ds)dr 
Jo 

■oo 

--< I H(t,s)f s (r)ds,h> r . 



We write this relation between H(t, r) and /t(r) in the form 

POO 

/ H(t,s)f s (r)ds = 5 t (r). 
Jo 



(47) 



Proposition 3.6 Let — 1 < 7 < 1, —1 < a < 1, and —1 < 7 + ct < 1. Then G 7 oG Q = G 1+a . 

Proof. The proof uses the following two relations: 

(1) G l9 (t) = J °°^ t H(t,s)g(s)ds, 7 e (-1,1); 

(2) J °° FyjHft, s)J r atS H(s,T)ds = T^f+ajH (t, r), with — 1 < 7, a < 1, and —1 < 7 + a < 1. 
Indeed, using (|46|) and changing the order of integration, we obtain the first relation 



G,g(t)= I ft(r)^(J H(r,s)g(s)ds)dr 

i roc 

g(s)( / f t (r)H(T,s)^dr)ds 



Tr r jH(t,s)g(s)ds. (48) 



It is readily seen that the internal integral in the second line of (148]) is meaningful, since 
/ t (r) is a function of exponential decay when r — > 00, which follows from (IT4"j) . Further, in 
order to show the second relation, we have 

poo poo poo poo 

/ ^ t H(t,s)^ s H(s,r)ds= (/ f t (p)H(p,s)p*dp)( f s (q)H(q,r)q a dq)ds 
Jo Jo Jo Jo 

POO 

ft(p)H(q,T)p\ a ( / H(p,s)f s (q)ds)dpdq. 



OO POO 



JO 



Due to (]4"T1) . this equals 

ft{p)p 1 {\ H( y q,r)q a 5 p ( y q)dq)dp = J H(p,r)p a f t (p)p^dp 

J 7 y+a,tH(t, 7~). 
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Now we are ready to prove the claimed semigroup property. Making use of the two proved 
relations, 



1. Proposition 13 . 6 1 immediately implies that G~ l = G_ 7 for arbitrary 7 G (—1,1). Indeed, 
G-y o G_ 7 = G = I, as well as G_ 7 o G 7 = /, where / is the identity operator. Thus, 
the statement in Corollary 13.51 is valid for all 7 G (— 1, 1). 

2. Proposition 13.61 remains valid for the family {Gt^, —1 < 7 < 1}, as well. 

3. The method used in this paper allows extension of results of Theorems l3.1H3.3l to the 
case of SDEs driven by time-changed linear fractional Levy stable motions. See [22] 
for CTRW limits of correlated random variables, whose limiting processes are time- 
changed fractional Brownian, or linear fractional Levy stable motions. 

4. The formula v(t,x) = Tu(t,x) for a solution of FPK equations associated with time- 
changed fBM, provides a useful tool for analysis of properties of a solution to initial 
value problems (l3~Tj) -( l32]) and (!37|) - (l38|) . as well as to the Cauchy problem for equation 



5. It is not necessary for the dependence of coefficients in (128]) on t to be of the form t 1 . 
This function can be replaced by [^(t)] 7 , where u{t) is a continuous function defined on 
[0,oo); however, the results essentially depend on the behavior of u{t) near zero and 
infinity. 
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